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ABSTRACT 



Aims. Among the known Ap stars, HD 101065 is probably one of the most interesting objects, demonstrating very rich spectra of 
rare-earth elements (REE). Strongly peculiar photometric parameters of this star that can not be fully reproduced by any modern 
theoretical calculations, even those accounting for realistic chemistry of its atmosphere. In this study we investigate a role of missing 
REE line opacity and construct a self-consistent atmospheric model based on accurate abundance and chemical stratification analysis. 
Methods. We employed the LLmodels stellar model atmosphere code together with DDAFit and Synthmag software packages to 
derive homogeneous and stratified abundances for 52 chemical elements and to construct a self-consistent model of HD 101065 
atmosphere. The opacity in REE lines is accounted for in details, by using up-to-date extensive theoretical calculations. 
Results. We show that REE elements play a key role in the radiative energy balance in the atmosphere of HD 101065, leading to 
the strong suppression of the Balmer jump and energy redistribution very different from that of normal stars. Introducing new line 
lists of REEs allowed us to reproduce, for the first time, spectral energy distribution of HD 101065 and achieve a better agreement 
between the unusually small observed Stromgren c\ index and the model predictions. Using combined photometric and spectroscopic 
approaches and based on the iterative procedure of abundance and stratification analysis we find effective temperature of HD 101065 
to be T cff = 6400 K. 
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1. Introduction 



Since its discovery by A. Przybylski in 1961 (Przybylski 19611) 
and named later after him, the Przybylski's star (HD 101065, 
hereafter PS) remains one of the most intriguing objects among 
Chemically Peculiar (CP) stars. These objects are usually char- 
acterized by abundance anomalies in their atmospheres and by 
non-uniform horizontal and vertical distributions of chemical el- 
ements and/or surface magnetic fields of different intensities. 

From many observed characteristics PS is a typical rep- 
resentative of the coolest part of rapidly oscillating A-type 
(roAp) stars. It pulsates with the typical roAp period, wl2 min 
dKurtz & Wegner|[l979l) . and it was the first roA p star discov- 
ered. It possesses a 2.3 kG surface magnetic field (ICowlev et al.l 
2000) and its low Fe abundance (an order of magnitude be- 
low the solar value) follows the trend of Fe abundance ver- 
sus effective temperature for Ap stars in 63 00-14000 K range 
(lKochukhovl2003HRvabchikova et ai1 l2004). A. Przybylski was 
the first to note a low Fe abundance a nd a strong overa bundance 
of some rare-earth (REE) elements dPrz vbvlskil 1 1 9661) . These 
features were co nfirmed in many subseque n t spectroscopic stud- 
ies of this star dWegner & Petfordll974l ICowlev et al.l 1 19771 
ICowlev & Mathvsl 119981: ICowlev et al.l 12000b . For instance, in 



Send offprint requests to: D. Shulyak, 
e-mail: denis . shulyak@gmail . com 

* Based on observations collected at the European Southern 
Observatory (Paranal, La Silla) and on data retrieved from the ESO 
Science Archive. 



the latter study ICowlev et all (l2000l) reported an overabundance 
by up to 4 dex (compared to the Sun) for elements heavier than 
Ni and enormously strong lines of the second ions of REEs com- 
pared to those of the first ions. This effect which is now com - 
monly known as REE anomaly (Ryabc hikova et all2001l 12004). 

As new data for REE transitions becomes available via labo- 
ratory measurements, more lines of REEs are identified and mea- 
sured in spectra of cool CP stars, thus improving abundance re- 
sults. Later on it it was discovered that PS is not even a champion 
in REE overa bundances. A t least a couple of hotter Ap stars: 
HP 17 0973 - iKatol d2003b ; HD 144897 - lRvabchikova"etal] 
(2006), with accurately derived REE abundances from numerous 
lines of the first and second ions have similar or even higher REE 
atmospheric overabundances. But a combination off a low tem- 
perature, low iron-peak abundances and large REE overabun- 
dances result in the extremely unusual observed spectrum of PS: 
weak lines of the iron-peak elements are lost in a forest of strong 
and numerous REE lines. For example, the typical density of 
only classified spectral features contributing to the theoretical 
spectra at 5375 A is 6-10 lines per A, and still not enough for 
proper description of the PS spectrum since every second feature 
can not be accounted for in the modern spectrum synthesis. 

Accurate spectroscopic analysis of PS requires dedicated 
model atmospheres. Indeed, abnormally strong REE absorption 
indicates that the atmospheric structure of the star may deviate 
significantly from that of normal, solar abundance stars. This is 
clearly seen, for instance, for the Stromgren c\ photometric in- 
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dex that is close to zero, which is very different than for normal 
stars. 

Previous attempt by iPiskunov & Kupkal d2001l) to account 
for peculiar absorption in the spectrum of PS and to fit color- 
indices based on the model atmosphere techniques faced a dif- 
ficulty of the absence of complete line lists of REE and the im- 
possibility to fit simultaneously photometric and spectroscopic 
features (like hydrogen lin es) with the same model atmos phere . 
IPiskunov & Kupk3d200Th used ATLAS9 (iKuruczll 19921 fl99l 
atmospheres with opacity distribution functions (ODF) recal- 
culated for individual abundances and artificially enhanced line 
(but not continuum) opacity of iron to simulate the missing REE 
opacity. They finally adopted T e ff = 6600 K, logg = 4.2 model 
as a preferable choice for spectroscopic analysis. However, this 
model failed to explain peculiar pho tometric parameters of the 
star. Enhanced metallicity model bv IPiskunov & Kup ka (2001) 
lead to the decrease of c\ index but its theoretical magnitude 
is still too high compared with the observed valu e. This model 
was u sed in extensive abundance study of PS by ICowlev et al.l 
(I2000h . 

Another striking evidence of the abnormality of the atmo- 
spheric structure of PS is the presence of the so-called core- win g 
anomaly in the hydrogen Balmer lines dCowlev et al.l 12001). 
This abrupt transition between the Doppler core of the hydro- 
gen lines and thei r Stark wings has been modeled empirically by 
Kochukh oTetail (120021) in terms of a temperature increase by 
500-1000 K at intermediate atmospheric heights but so far has 
not been explained theoretically. 

Thus, HD 101065 remains an extremely challenging and 
very intriguing target for the application of advanced model 
atmosphere analysis and high-resolution spectroscopy. In this 
paper we present a new attempt to combine both of these ap- 
proaches and to construct a self-consistent model atmosphere of 
PS taking into account realistic chemistry of its atmosphere and 
employing up-to-date theoretical calculations of the spectra of 
selected REEs. 

Starting with a brief description of observations in Sect. [2] 
we then present methods of analysis in Sect. [3] and describe the 
new lists of REE lines employed in the model atmosphere calcu- 
lations. Results of our analysis are presented in Sect. [4] and are 
followed by general conclusions in Sect. [5] Our paper concludes 
with a brief discussion in Sect. [6] 



3. Methods 

3. 1 . Abundance analysis 

The main goal of the present paper is to calculate a model at- 
mosphere of PS accounting for its anomalous chemical compo- 
sition. In magnetic chemically peculiar star an accurate abun- 
dance analysis may be carried out by fitting the magnetic syn- 
thetic spectrum to the observed line profiles. This procedure is 
very time consuming. Therefore, to simplify calculations, we 
analysed the equivalent widths replacing magnetic effects by 
a pseudo-microturbulence. It is sufficient for line opacity cal- 
culations, although the abundances derived using the common 
value of microturbu lent velocity m ay be slightly inaccurate (see 
discu ssion i n ICowlev et alj|200 0). In present analysis we used 
fhelKuruczl (Il993l) WIDT H9 code modified by V. Tsymbal (see 
iRvabchikova et allll997h . 

The VALD database (IPiskunov et all Il995t iKupka et al l 

1999) and the DREAM REE line database dBiemont et all 1999. 
and references therein), which is made accessible via the VALD 
extraction procedures, are our main sources for atomic param- 
eters. For the REE, which provide the majority of lines in the 
PS spectrum, accurate laboratory tran sition probabi l ities m ea- 
sured by Wiscons i n grou p are u sed: lLawler etal] (|2001al) - 
Lan; lLawler et al.l (l2001al) - La n : lLawler et alj d2009b - Ce n: 
| Den Hartog et al.l (120031) - N dn; lLawler et alj d2006l) - Sinn; 
Lawler et al. l d2001bl) - Eu n: iDen Hartog et all d2006l) - Gdn; 
Lawler et al d2001d) - Tb n ; IWickliffe et all d2000t) -~5 v i/Dy n; 
Lawler et al (p008) - Er n; IWickliffe & Lawleri d 19971) - Tm n; 
lLawler et all d2007l) - Hfn. Among heavier elements Th and 
U are particularly interesting. For these species we used lab- 
oratory and calculated transition probabilit ies from the f ollow- 

Thn: iBiemont et all d2002) 
Un. Transition probabilities 



Nilssonetal. (2002a) 



ing papers: 

- Thin; iNilsson et all (I2002bl) 



for REEs in the second ioniz ation stage were taken fr om the 
DREAM database except Pr m iMashonkin a et all (120091) . Nd m 
IRvabchikova et all d2006l) . Eu m IWvart et all (l2008l) 7Tb m and 
Dy in (Ryabtsev, private communication). 

Extended measurements and theoretical calculations for the 
REEs, Th and U justify revision of some of the partition func- 
tions (PF), which until now mainly came from iKura cz (1993) 
ATLAS9 code. Our recalculated PFs are made available onlineQ. 

We did not take into account the hyperfine structure (hfs) in 
our abundance analysis because the hfs constants are available 
only for a small subset of spectral lines. 



2. Observations 

In our study we employed several existing spectra of PS: NTT 
spectrum with the resolving power R «s 80 000 in the 3950 - 
6630 A wavelength region, dCowlev et al .1120001) . UVES spec- 
trum in the 3 100-10 000 A region with the resolving powe r R * 
80000 dKochukhov et all 120021: IRvabchikova et alJl2008l) . and 
an averaged UVES spectrum from the time-series observations 
obtained in the 4960 -6990 A spectral region with the resolving 
powe r R « 1 15 000 dKochukhov et al.ll2007t IRvabchikova et all 
[2007h . 

Details of the spectroscopic observations and full descrip- 
tion of the data reduction are given in the corresponding papers. 
Direct comparison of the normalized spectra reveals a remark- 
able agreement between all observations. This allows u s to us e 
the equivalent width measurements from ICowlev et alJ (yOOO), 
supplementing them with the measurements in the spectral re- 
gion A > 6600 A using the other two spectra. 



3.2. Stratification analysis 

In the majority of cool Ap stars the atmospheres are not chem- 
ically homogeneous. Element stratification built up by atomic 
diffusion influences the atmospheric structure dShulvak et al] 
2009; Kochuk ho7etalir2 009: Leb lanc et al]l2009l) . As it follows 
from Kh an & Shuly ak (2007), Fe, Si, and Cr are the elements 
that have the major effect on the atmospheric T - P structure in 
the case of homogenous abundances. The cumulative effect of el- 
ement stratificatio n on the mod el structure was demostrat ed by 
IWade etall d2003l) and later by iMonin & LeBland d2007l) who 
used diffusion calculations for 39 chemical elements simultane- 
ously. Thus, the stratification analysis of PS was performed for 
four elements: Fe, Ba, Ca, and Si. It was not possible to carry 
out stratification analysis for Cr and other iron-peak elements 
due to weakness of their absorption features and heavy blend- 
ing by REE lines. The lower energy levels of practically all Cr 

1 http : //www . astro . uu . se/~ oleg/pf . html 
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lines used in the abundance determination belong to a narrow 
energy range, which makes the Cr lines insensitive to abundance 
gradients. Hence, we restricted stratification analysis to the four 
elements since they are represented in HD 101065 by a sufficient 
number of atomic lines, probing different atmospheric layers and 
enabling line profile fitting. 

We applied a step-function approximation of chemical strat- 
ification as implemented in DDAFit - an automatic procedure 
for de termination of vertical abundance gradients ( Koch ukhovl 
120071) . In this routine, the vertical abundance distribution of an 
element is described by four parameters: chemical abundance in 
the upper atmosphere, abundance in deep layers, the position of 
the abundance jump and its width. All four parameters can be op- 
timized simultaneously with the least-squares fitting procedure. 

Most unblended lines in the PS spectrum are located in the 
red spectral region, where line profiles are also substantially dis- 
torted by the Zeeman splitting due to wavelength dependence 
of the Zeeman effect. DDAFit enables an accurate stratification 
analysis of such magnetically-splitted lines. In the search for op- 
timal vertical distribution of chemical eleme nts we use magneti c 
spectrum synthesis with the Synthmag code (Kochukhov 2007). 
This softw are represents an improved version of the program de- 
veloped bv lPiskunovl (1 19991) . 

A list of the lines used in stratification analysis is given in 
Table [3] (Online material). To get a better sensitivity of the fit- 
ted Ca distribution to the upper atmospheric layers, we used the 
IR-trip let Ca n line at A 8498 A. As shown by Rvabchi kova et al.l 
(2008J) and lCowlev et al.l d2009l) . its core is represented entirely 
by the absorption of the heavy isotope 48 Ca. Therefore wave- 
length of the line of this isotope is given in Table [3] 



3.3. Calculation of model atmospheres 

To perform the model atmos phere calculations w e used the re- 
cent version of the LLmodels ( Shulvak et al. 2004) stellar model 
atmosphere code. The code accounts for the effects of individ- 
ual and stratified abundances. The stratification of chemical ele- 
ments is an input parameter and thus does not change during the 
model atmosphere calculation process. This allows us to explore 
the changes in the model structure due to stratification that was 
inferred directly from observations, without modeling poorly un- 
derstood processes that could be responsible for the observed 
inhomogeneities. 

Note that such an empirical analysis of chemical element 
stratification is based on the model atmosphere technique and 
thus the temperature-pressure structure of the model atmosphere 
itself depends upon the stratification which we want to deter- 
mine. Therefore, the calculation of the model atmosphere and 
the stratification (abundances) analysis are linked together and 
an iterative procedure should be used in this case. It consists 
of repeated steps of stratification and abundance analysis that 
provide an input for the calculations of model atmosphere until 
atmospheric parameters (T e ff, logg, etc.) converge and theoreti- 
cal observables (photometric colors, profiles of hydrogen lines, 
etc.) fit observations. We refer the reader t o the recent papers 
dShulvak et al.l2009tlkochukhov et al.l2009t) for a more detailed 
description of this technique. 

The VALD database dPiskunov et all 119951: iKupka et al.l 

1 19991) was used as a main source of the atomic line data for 
computation of the line absorption coefficient. The recent VALD 
compilation contains information on about 66 x 10 6 atomic tran- 



sitions. Most of them come from the latest theoretical calcula- 
tions performed by R. KurictQ. 

Since a 2.3 kG surface magnetic field of PS is too weak 
to affect no ticeably the atmospheric structure and energy d istri- 
bution (see iKochukhov et allEooH iKhan & Shulvakll2006L for 
more details) we adopted a 1 km s _1 micro turbulent velocity to 
roughly account for the magnetic broadening and to avoid time- 
consuming model calculations with detailed polarized radiative 
transfer. 



3.4. REE line lists for model atmosphere calculation 

The main sources for the REEs transition probabilities data 
between experimentally known energy levels are VALD and 
DREAM databases, hereafter referred to "VALD data". These 
data were supplemented by the transition pr obabilities for the 
obser ved and pred icted lines of Prn-Prm ( Mashonkina et al 
2009), Nd n- Nd m dMashonkina et al.ll2005l iRvabchikova et al 



2006), Eu in dWvart et al.ll2008l) . Tb m and Dy in (Ryabtsev, pri- 
vate communication), hereafter referred to as "ISAN data". For 
the latter two ions transition probability calculations were based 
on the extended term analysis. 

Table|2]demonstrates a dramatic difference between the num- 
ber of observed (VALD) and predicted (ISAN) transitions for Pr 
and Nd, in particular for the lines of the first ions. We are miss- 
ing more than 90% of the potential line absorbers. The same 
situation is expected for lines of the other REEs. Taking this 
situation into account, we performed calculations for Smn - 
the second-abundant element after Nd in PS. The energy lev- 
els of Sm ii we re calcu l ated b y the Hartree-Fock method imple- 
mented in the ICowanl d!981l) code. The ground state of Smn 
is the 4f*6s configuration. In addition, our calculations include 
even configurations 4f 6 nd (n=5,6), 4^78, 4f 5 5d6p, 4f 5 6s6p, and 
odd configurations 4f 6 np (n=6-8), 4f 4 5d 2 6p, 4^68^, 4f 5 5d 2 , 
4f 5 5d6s. Calculations are based on the wave functions obtained 
by the fittings the energy levels. As in the case of Pr n and Ndn, 
all Hartree-Fock transition integrals are scaled by a factor 0.85. 
In total, transition probabilities for more than one million lines 
were calculated and added to the ISAN list of the REE l in es. 

The recent studies by iMashonkina et al.l (120051 120091) 
demonstrated that the line formation of Pr and Nd can strongly 
deviate from the local thermodynamic equilibrium. For instance, 
the doubly ionized lines of these elements are unusually strong 
due to combined effects of stratification of these elements and 
departures from LTE. However, a precise NLTE analysis is be- 
yond the scope of this paper and currently could not be coupled 
to a mod el calculation. Theref ore, we followed the approach out- 
lined in IShulvak et al.l (120091) where authors used a simplified 
treatment of the REE NLTE opacity. 

Taking into account a systematic difference in abundances 
derived for the first and second REE ions, it is essential to re- 
duce the oscillator strengths for the singly ionized REE lines 
by this difference in the model line list while using the abun- 
dances derived from second ions as an input for model atmo- 
sphere calculations (or vice versa). Obviously, the adopted re- 
duction factors depend upon the assumed abundances of the first 
and second ions and thus change slightly if the model param- 
eters (r e (j, logg, etc.) are modified in the course of the itera- 
tive procedure of abundance analysis. The logarithmic scaling 
factors, log(A^n/A^in), adopted for our final model are (in dex): 
-1.6 (Ce), -2.4 (Pr), -1 (Nd), -2.66 (Tb), -2.5 (Dy), and -1.64 
(Eu). These factors were applied to all lines of the respective 
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ions in the master line list. Then this line list was processed by 
the line pre-selection procedure in the LLmodels code to select 
only those lines that contribute non-negligibly to the total line 
opacity coefficient for a given temperature-pressure distribution. 
Applying this line scaling procedure allowed us to mimic the line 
strengths correspond to the NLTE ionizatio n equilibrium abun- 
dance for each REE (see lShulvak et al.|l2009l for more details). 



4. Results 

4. 1 . Adopted atmospheric parameters 

We started analysis with a homogeneous abundance model cal- 
culated with T e ff = 6600 K, logg = 4 -2, and individual abun- 
dances taken from lCowlev et all ( 2000I) . Once the stratification 
was introduced in the model atmosphere, the fit to the hydro- 
gen Ha line and the observed photometric parameters required 
decreasing of T e ff of the star down to approximately 6400 K. In 
total, four iterations of the abundance analysis and the model 
atmosphere calculation were performed to achieve a converged 
solution. The abundances from the third iteration did not intro- 
duce noticeable changes in the atmospheric structure. 

The effective temperatu re obtained in our stud y is identi- 
cal to the T e ff obtained by iKochukhov et all (120021) . who fitted 
the Ha and H/3 lines using the same UVES spectra. However, 
Kochukhov et al. relied on less advanced ATLAS9 model atmo- 
spheres cal culated using ODFs with a rtificially increased Fe line 
absorption (Piskunov & Kupka 2001) and also attempted to em- 
pirically adjust the T-t structure of their model. 

Figure Q] illustrates the observed and predicted Ha 
line profiles calculated with different atmospheric models. 
Unfortunately, it is impossible to infer an accurate value of log g 
from the Ha profile due to low temperature of the star and re- 
sulting poor sensitivity of the hydrogen lines to the pressure 
stratification. For this reason we kept this parameter fixed dur- 
ing the iterative procedure of abundance analysis (however, see 
Discussion). 



1.2 



4.2. Abundance pattern 

Table |4] (Online material) summarizes individual abundances of 
chemical elements in the atmosphere of PS derived with the final 
model atmosphere. Note the strong enrichment of REE elements 
compared to the solar composition. These elements dominate the 
spectrum of PS, playing a key role in the radiative energy bal- 
ance (see below). 

In Fig. [2] we illustrate abundance pattern of PS relative to 
the re cent compilation of the solar abundances dGrevesse et al.l 
l2007h . The elements heavier than Ba exhibit an overabundance 
by 3 dex and larger. Several REEs show clear ionization anomaly 
related to stratification of these elements and departures from 
LTE. 

Vertical distribution of Si, Fe, Ba, and Ca is presented in 
Fig. [3] This figure illustrates the difference between the initial 
and final stratification as derived at the first and the last iterations 
respectively. Ba seems to be the only element whose stratifica- 
tion did not change appreciably during iterations. Only the posi- 
tion of the abundance jump changed slightly. Fe shows a change 
in the abundance in the upper atmosphere by about ~ 1 dex. The 
position of the abundance jump for Si changed significantly. In 
contrast, the Ca abundance in the upper atmosphere changes dra- 
matically. These results illustrate the importance of taking into 
account a feedback of stratification on the model structure. 
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Fig. 1. Comparison between the observed and predicted Ha line 
profiles. 
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Fig. 2. Chemical composition of the atmosphere of HD 101065. 
Different symbols correspond to ions from I to EH. Bars show 
abundance ranges for stratified elements. 



The effect of introducing stratification on the fit to Ha line 
appears to be not very strong. The maximum difference between 
the stratified and non-stratified model predictions is at the level 
of 1% as seen from Fig. Q] However, the two profiles are still 
clearly distinguishable. Thus, chemical stratification should be 
taken into account in accurate modeling of the hydrogen line 
formation. 

The Fe and Ba abundance distributions agree with those de- 
rived previously by either trial-and-error met hod ([Rvabchikova 
2003) or based on the equivalent widths (lYushchenko et al. 
2007). Figures |3]|6] (Online material) illustrate the resulting fit 
to line profiles. 
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Fig. 7. Observed and theoretical energy distribution of HD 101065 calculated using line lists from two datasets: VALD and ISAN. 
For comparison purpose the theoretical energy distribution of a normal F-type star HD 49933 is also shown. All theoretical calcula- 
tions are smoothed with 20A Gaussian. 
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Fig. 3. Stratification of four elements in the atmosphere of 
HD 101065 derived at the first and last iteration of the abun- 
dance analysis. 



4.3. Energy distribution and photometric colors 

Unusual photometric parameters of PS clearly indicate that its 
energy distribution should differ much from that of normal stars 
with similar r e ff . Indeed, a dense forest of REE lines obtained 
from theoretical computations should have a strong impact on 
the atmospheric energy balance. 

Figure Q illustrates theoretical energy distributions of 
HD 101065 and a normal F-type star HD 49933. Both stars 
have close effective temperatures but very different atmospheric 
chemistry: HD 49933 has a slight underabundance of iron-peak 
elements and a ne a rly so lar abundance of REEs as derived by 
Rvabchikova et al. (2009). Model fluxes were computed with 
the stratified distribution of Si, Ca, Fe, and Ba as well as with 
the scaled REE abundan ces as described above. I n addition, 
ATLAS9 ODF model from lPiskunov & Kupkal(l200lh calculated 
with T e ff = 6600 K, logg = 4.2 and Fe opacity scaled to simu- 
late an enhanced REE absorption is shown for comparison pur- 



pose. To keep the figure representative, here we show only mod- 
els with final T e ff = 6400 K, log g = 4.2 and different assump- 
tions about REE opacity. Energy distributions of models with 
r eff = 6500 K and T eff = 6600 K are displayed in Fig.[8](Online 
Material). 

The observed energy distribution of PS is constructed com- 
bining the large-aperture IUE observations extracted from the 
INES database^ and pho tometric data i n the UBV (|Wegner 
119761). Genevffl JH K dCatalano et al l 1 19981) and 2MASS 
(ISkrutskie et al.ll2006l) systems converted to absolute fluxes. The 
model fluxes are scaled to account for the distance to PS, which 
is 112 + 11 pc accordi ng to the revised H ipparcos parallax of 
7r = 8.93 + 0.87 mas (Ivan Leeuwenl 12007). Theoretical fluxes 
are further scaled to match the observed Paschen continuum de- 
fined by the UBV and Geneva photometry. This scaling fac- 
tor corresponds to the radius of the star, for which we found 
R — 1 .98 + 0.03/? o . The same scaling is applied to all theoretical 
fluxes presented in Fig. [7] For models with T e ff - 6500 K and 
r eff = 6600 K we find R = 1.9O/? and R = 1.82tf respectively 
(see Fig. [8] Online Material). 

It is evident that the enormous REE absorption leads to 
smoothing of the Balmer jump. This is clearly visible for both 
models computed with the new ISAN REE line list and with the 
original VALD REE line list. The only difference is that the for- 
mer one predicts less flux redward of 4 100 A and more flux in the 
region 3000A and 3800A (note that the model computed with 
mean REE abundances should be used with big caution due to 
ignoring REE anomaly, and we present it in Fig.|7]for illustrative 
purpose only). Taking into account difficulties in predicting ac- 
curately the spectra of REEs, one can conclude that our models 
computed with the recent improvements of the line REE opac- 
ity provide a reasonably good fit to the observed fluxes. They 
are able to reproduce not only the shape of the observed energy 
distribution but also the amplitude of the Balmer jump reason- 
ably well, much better than in the previo us attempt by ODF- 
based modeling dPiskunov & Kupkal2001l) . However, there are 
still some discrepancies with observations, i.e. too high fluxes in 

3 http : // sdc . laef f . inta . es/ ines/index2 . html 

4 ht tp : //ob swww . uni g e . ch/ g cpd/ph 1 3 . html 
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the 2000-2500 A region and around 4500 A. The description of 
the energy distribution of PS can be probably improved in the 
future once more complete line lists of REEs will become avail- 
able. A new well-calibrated observed spectrophotometry in the 
3300-10000 A region would also be extremely useful for further 
modeling. 

The impact of the peculiar energy distribution on the pho- 
tometric parameters of HD 101065 is illustrated in Table Q] 
Observations in the Stro mgren system were taken from 
Hauck & Mermillioi (fT998h . Predictions of Pisk unov & K upka 
d2001l) model calculated wit h the PDF line opac ity representa- 
tion using abundances from dCowlev et alJ feOOO) are also given 
for comparison. We also show theoretical photometric parame- 
ters for the models calculated with the mean REE abundances 
corresponding to the first and second ions, as well as for the 
models with the scaled REE opacity, and for different T e g. All 
models except the ODF-based one were computed with the 
chemical stratification illustrated in Fig. [3] 

There are few things to note from Table Q] First of all, a 
strong impact of the REE opacity is clearly seen for the C[ photo- 
metric index which reaches an even more negative value than the 
observed one for the models with the ISAN data and the mean 
REE abundances. Similarly, the recent version of VALD already 
provides REE opacity that is enough to bring down ci index sig- 
nificantly compared to previous calculations and provides a re- 
sult very close to the observed one. However, mean abundances 
can not be used for accurate modeling due to significant REE 
ionization anomaly. Accounting for this via the scaling of REE 
line opacity decreases the flux redistribution effect and thus the 
impact on the c\ index. The implementation of the new ISAN 
data gives a much better fit to observations than using the data 
from VALD alone, yet this fit is far from being fully satisfactory. 

Second, it remains difficult to obtain a good fit simultane- 
ously for all photometric indices. For instance, all models fail 
to reproduce the U - B value, while b - y is not reproduced by 
models with the scaled REE opacity unless one decreases T e g 
to as 6200 K. However, lowering r e jf is not consistent with the 
observed Ha line. Similar, worse fit is obtained for models with 
T e s > 6400 K, which show reduced b — y and increased ci indi- 
cies. 

Third, stratification of Fe, Si, Ca, and Ba considered in our 
study does not play a significant role for the overall energy re- 
distribution and thus photometric colors as seen from the sixth 
raw of Table Q] Most of the photometric indicators are negligibly 
affected by the presence of stratification. The c\ index is again 
an exception since it decreases by * 0.03 mag if stratification is 
introduced. However, this change is still fairly small compared 
to typical precision of the observed Stromgren photometry. 

Finally, a high sensitivity of the photometric parameters to 
the new ISAN calculations indicates that the general disagree- 
ment between models and observations is the subject of future 
more precise and complete calculations of REE spectra. Using 
the data for more REE elements than those presented in this in- 
vestigation can potentially improve the fit to the highly peculiar 
observed photometric properties of PS. Implementation of the 
ISAN line lists for only a few REE ions already allowed us to 
achieve significantly better results compared to previous model- 
ing attempts and clearly showed that the "enigmatic" character- 
istics of PS likely result from the REE line opacity which was 
ignored in previous calculations. We can thus suggest that fur- 
ther step in the modeling of the observed properties of PS should 
be made in the direction of improved theoretical calculations of 
REE transitions. 



Table 1. Observed and predicted colors of HD 101065. 





b - y 


mi 




B - V 


U - B 


Observations 


0.442 


0.427 


-0.013 


0.760 


0.200 


t6600g4.2, ODF 


0.387 


0.582 


0.298 


0.767 


0.571 


t6400g4.2 
VALD 

REE mean abun. 


0.440 


0.439 


0.029 


0.745 


0.268 


t6400g4.2 
ISAN 

REE mean abun. 


0.456 


0.587 


-0.202 


0.829 


0.379 


t6400g4.2 
VALD 

REE scaled abun. 


0.394 


0.420 


0.312 


0.681 


0.368 


t6400g4.2 
ISAN 

REE scaled abun. 


0.384 


0.413 


0.158 


0.663 


0.274 


t6500g4.2 
ISAN 

REE scaled abun. 


0.367 


0.397 


0.200 


0.634 


0.255 


t6600g4.2 
ISAN 

REE scaled abun. 


0.350 


0.391 


0.229 


0.607 


0.240 


t6400g4.2 
ISAN, horn. 
REE scaled abun. 


0.383 


0.396 


0.188 


0.656 


0.271 


t6200g4.5 
ISAN 

REE scaled abun. 


0.417 


0.444 


0.024 


0.719 


0.261 



All models were computed with stratification of Fe, Si, Ca, and Ba, 
except those marked as "ODF" and "horn." 



The strong suppression of the Balmer jump and energy re- 
distribution to the infrared is a direct result of the unusually high 
REE abundances in the stellar atmosphere. Table [2] summarizes 
some statistics illustrating importance of the REE opacity in the 
atmosphere of PS. It shows the number of lines of a given ion 
before and after line preselection procedure performed by the 
LLmodels code. This procedure selects from the master line list 
only those lines that noticeably contribute to the line opacity 
for a given T-P structure. Usually the preselection criterion is 
Kline /tfcont > 1%, where Kii ne and K cont are the line center and 
continuum opacity coefficients, respectively. It is seen that, for 
instance, all lines of Pr n and Nd n present in the current version 
of VALD contribute to the opacity in the atmosphere of PS. A 
larger number of lines is selected from the more complete ISAN 
line lists. 

As expected from the low magnitude of the c\ index and the 
energy distribution shown in Fig. [7] a great majority of strong 
REE lines are concentrated around the Balmer jump. Figure [9] 
illustrates the line distribution for selected REE ions depending 
upon their wavelength position and central intensity. We show 
the number of lines per 50 A wavelength bin for which the cen- 
tral intensities are greater or equal to a given normalized inten- 
sity. The latter ranges from (fully saturated strong line) to 1 
(continuum level, weak line). One can seen that ions like Cen, 
Prn, Ndn/m, Smu have a large number of strong lines located 
exactly in the region of Balmer jump and blueward, dramatically 
contributing to the opacity and radiation field. 
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Table 2. Summary of REE lines in VALD and combined ISAN line lists before and after line preselection procedure. 



Pru 



Prm 



Nd ii Nd m 



Smn 



Eu in Dy m Ce n 



Cem 



VALD 
ISAN 



510 

103428 



19104 



1285 
1172579 



Original 
71 1346 
6858 1047088 



893 
23827 



1313 
31121 



16046 2966 



After preselection with 1% limit for T tl 



VALD 
ISAN 



510 
83167 



6539 



1285 
74413 



= 6400 K, logg = 4.2 model, REE scaled abundances 

71 1121 58 1134 15908 1270 
3810 58637 74 21715 



N 




Wavelength, (nm) Wavelength, (nm) Wavelength, (nm) 



N 




Wavelength, (nm) Wavelength, (nm) Wavelength, (nm) 



Fig. 9. Distribution of lines for several REE ions depending on their position and central depth in the ISAN line lists for the model 
atmosphere with T e ff = 6400 K, logg = 4.2. 



5. Conclusions 

HD 101065 is known to be one of the most peculiar stars ever 
found. Its unusual photometric colors and extremely rich spectra 
of REE elements represent a major challenge for modern stellar 
atmosphere theory. With the increasing number of the accurately 
measured and calculated REE transitions and the progress in the 
model atmosphere techniques it became possible to carry out this 
quantitatively new analysis of the spectrum of HD 101065 and 
to re-determine atmospheric parameters of this star. 

In this study we investigated the impact of the new theo- 
retically computed spectra of some REE elements on the atmo- 
spheric properties of HD 101065 and carried out a detailed spec- 
troscopic analysis of the stratification of Si, Fe, Ba, and Ca in 
the stellar atmosphere. Theoretical computations of REE transi- 
tions have extended already existing line lists of REE lines by a 
factor of one thousand. Implementing an iterative procedure of 
the abundance analysis, we re-derived atmospheric parameters 
of the star in a self-consistent way, accounting for the ionization 
disequilibrium between the first and second ions of the REEs 
caused by significant NLTE effects. 

The main results of our investigation can be summarized in 
the following: 

- REEs appear to play a key role in the radiative energy bal- 
ance in the atmosphere of HD 101065, leading to the strong 



suppression of the Balmer jump and energy redistribution 
from UV and visual to IR. 

- The abnormal photometric color-indices of HD 101065 arise 
due to a strong absorption in REE lines. We showed that 
the introduction of the new extensive line lists of REEs al- 
low one to achieve a better agreement between the unusually 
low magnitude of the Stromgren c\ index observed for this 
star and the model predictions. We also demonstrated, for the 
first time, a satisfactory agreement between the observed and 
computed spectral energy distribution of HD 101065. 

- The remaining discrepancy between the observed and cal- 
culated photometric parameters indicates that some impor- 
tant opacity sources are still missing in the modeling process. 
For this reason the future extension of the REE line lists is 
needed for complete understanding of the observed proper- 
ties of HD 101065. 

- Using theoretical spectrum synthesis we derived stratifica- 
tion profiles for Si, Ca, Fe, and Ba. We found a strong de- 
pletion of all these elements in the upper atmosphere of 
HD 101065. The effect of stratification on the Ha line profile 
appears to be at the level of 1 %, which is however important 
for accurate calculation of the hydrogen line profiles. 

- There is not much influence of the stratification on the con- 
sidered photometric parameters of HD 101065, except for 
the Stromregn c\ index, which is modified by stratification 
by 0.03 mag. 
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- Using combined photometric and spectroscopic analysis and 
based on our novel iterative procedure of the abundance and 
stratification analysis we find effective temperature of the 
HD 1 1 065 to be r eff = 6400 K. 

Taking into account the current progress in the atmospheric 
modeling, understanding of the REE opacity, and detailed spec- 
troscopic analysis we conclude that there is nothing special in 
HD 101065 compared to other known cool Ap and roAp stars 
except its relatively low r e ff and a high abundance of REEs. 
Consequently, for this moment HD 101065 appears to be the 
coolest Ap and roAp star known but otherwise is a typical mem- 
ber of these groups. 



6. Discussion 

We can compare the global parameters of HD 101065 derived 
in the present study with the values obtaine d by an independent 
pulsation modeling ( Mk rtichian et alj |2008). The best model ex- 
plaining the observed frequency pattern of HD 101065 is char- 
acterized by r eff = 6622+100 K, log g = 4.06+0.04, R= 1.90 + 
0.08 R e . These parameters agree with ours within 2<x. However, 
to fit the observed frequencies in the framework of current as- 
terosismic models one needs dipolar magnetic field strength 
Bp = 8.7 kG, which is incompatible with the observed field mod- 
ulus of (B) = 2.3 kG. One of the reasons for this significant dis- 
crepancy between the observed and inferred magnetic field is 
the use of normal T-t relation in the pulsation modeling. This 
relation differs significantly from a more realistic one derived in 
the present paper, as illustrated in Fig. [10] (Note that comput- 
ing model with solar abundances we did not apply scaling of the 
REE opacity due to the low abundance of REE elements in solar 
atmosphere and thus its negligible impact on model structure). 
As expected from the fit to Ha line and photometric parameters 
presented above, the effect of stratification on model T — P struc- 
ture is negligible and hardly distinguishable from the homoge- 
neous abundance model (dash-dotted line in Fig.fTOt. Note that it 
is still considereably different from the solar abundance model, 
and is in good agreement with the result of M onin & LeB lanc 
(2007). 

In spite of our new results, illustrating an improved under- 
standing of the physics of the atmosphere of HD 101065, it is 
still impossible to describe all observed properties of this star 
with the current state-of-the-art modeling. In this section we 
would like to discuss physical effects that we ignored in our anal- 
ysis. 

Taking into account a jungle of the REE lines seen in the 
spectrum of HD 101065 and their importance as an opacity 
source for the model atmosphere calculation it would be rea- 
sonable to search for the possible stratification of rare-earth el- 
ements. Considering the case of a hotter roAp star HD 24712, 
for which the accumulation of REEs in the upper atmosphere 
plays a critical role , producing a charac teristic inverse temper- 
ature gradient (see IShulvak et al.l 120091 for more details), one 
would expect the same mechanism to operate in the atmosphere 
of HD 101065 as well. However, the stratification analysis of 
REEs is highly complicated due to blending of REE lines with 
accurately known atomic parameters (that ideally should be used 
for such an analysis) by other REE lines for that no laboratory 
measurements exist and thus can not be included in the spectrum 
synthesis. At the same time, assuming that REEs are distributed 
similarly in both stars, we do not expect a strong impact of the 
REE stratification on the energy distribution of HD 101065 sim- 
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Fig. 10. Comparison in atmospheric temperature structure 
of theoretical mo d els. T hick line - ODF model from 
iPiskunov & Kupkal (l200ll) . dotted line - homogeneos solar 
abundances model, dashed line - stratified abundances final 
model adopted in the present study, dash-dotted line - same as 
before but assuming homogeneos abundances. 



ply because it would affect only the optically thin atmospheric 
layers. 

As mentioned above, due to the low T e ff of PS, it is impossi- 
ble to infer the log g value from the hydrogen lines. Furthermore, 
such logg indicator as the c\ photometric index is strongly af- 
fected by the REE opacity, which is still far from being com- 
pletely understood. The ISAN calculations gave us information 
on several REE ions but the role of other REE species is not 
known at this moment. On the other hand, as seen from the high 
resolution spectra of HD 101065, most of the line absorption is 
due to such ions as Prii/m, Ndii/m, and Ceii/m that are repre- 
sented in our calculations by a large number of lines. This sug- 
gests that the magnitude of the energy redistribution is well mod- 
eled. However, extending our calculations using new line lists for 
other REE ions will allow to improve the accuracy of theoretical 
c\ index and constrain the \ogg of the star. In principle, increas- 
ing log g up to 4.5 and decreasing T e $ down to 6200 K provides 
a better fit to the observed photometric parameters as demon- 
strated in the last row of Table \T\ Also the fit to the Ha line 
remains reasonable with these parameters. However, the ioniza- 
tion equilibrium for different elements is not as good as with the 
T e ff = 6400 K, logg = 4.2 model. This is why, and taking into 
account various uncertainties related to the REE opacity, we did 
not attempt to derive surface gravity in our study. This limitation 
does not affect any of the results presented in this paper. 
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Wavelength 

Fig. 3. Comparison between the observed and theoretical profiles 
for the Ba lines. Observations are shown with symbols, thick 
full line corresponds to the synthetic spectrum calculated with 
the best-fitting stratified abundance distribution and dashed line 
shows prediction of the spectrum synthesis with a homogeneous 
abundance. 
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Fig. 4. Same as in Fig[3]but for Ca. 



Fig. 5. Same as in Fig[3]but for Si. 
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Fig. 6. Same as in Fig[3]but for Fe. 
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Table 3. List of Si, Ca, Fe, and Ba spectral lines used for reconstruction of chemical stratification. 



Ion 


A (A) 


Ei (eV) 


log(g/) 


log 7s, 


Ref. 


Si i 


5517.533 


5.082 


-2.49 


-4.46 


solar 


Si i 


5690.425 


4.930 


-1.76 


-4.57 


NIST 


Si i 


5797.688 


5.614 


-2.67 


-3.86 


K07 


Si i 


5797.856 


4.954 


-2.05 


-4.32 


NIST 


Si i 


6087.697 


5.871 


-2.38 


approx 


K07 


Si i 


6087.805 


5.871 


-1.71 


approx 


solar 


Si i 


6155.134 


5.619 


-0.80 


-3.70 


solar 


Sin 


6347.109 


8.121 


0.297 


-5.04 


BBCB 


Sin 


6371.371 


8.121 


-0.003 


-5.04 


BBCB 


Si i 


8742.446 


5.871 


-0.510 


-4.52 


solar 


Si i 


8892.720 


5.984 


-0.830 


-4.37 


solar 


Cai 


6717.681 


2.709 


-0.524 


-4.90 


SR 


Cai 


7148.150 


2.709 


0.137 


-6.01 


SR 


Cai 


7326.145 


2.933 


-0.208 


-5.16 


S 


Can 


8254.721 


7.515 


-0.398 


-4.60 


T 


Can 


8498.223 


1.692 


-1.416 


-5.70 


T 


Can 


8912.068 


7.047 


0.637 


-5.10 


T 


Fei 


4404.750 


1.557 


-0.142 


-6.20 


MFW 


Fen 


5169.033 


2.891 


-1.14 


-6.59 


M 


Fei 


5405.775 


0.990 


-1.844 


-6.30 


MFW 


Fei 


5415.199 


4.386 


0.642 


-4.76 


BWL 


Fei 


5615.644 


3.332 


0.050 


-5.50 


BKK 


Fei 


6219.281 


2.198 


-2.433 


-6.20 


MFW 


Fen 


7222.394 


3.889 


-3.36 


-6.67 


HLGH 


Fei 


8248.129 


4.371 


-0.887 


-5.36 


K07 


Fei 


8365.634 


3.251 


-2.047 


-6.10 


BK 


Fei 


8945.189 


5.033 


-0.235 


-5.17 


K07 


Bai 


5777.618 


1.676 


0.72 


approx 


MW 


Ban 


5853.668 


0.604 


-1.00 


-6.35 


MW 


Bai 


6019.465 


1.120 


-0.10 


approx 


MW 


Bai 


6063.109 


1.143 


0.21 


approx 


MW 


Ban 


6141.713 


0.704 


-0.076 


-6.35 


MW 


Ban 


6496.897 


0.604 


-0.377 


-6.35 


MW 


Bai 


7059.938 


1.190 


0.68 


approx 


MW 


Bai 


7280.297 


1.143 


0.47 


approx 


MW 



Columns give the ion identification, central wavelength, excitation potential, oscillator strength (log(g/)) and the Stark damping constant (log yst) 

for T = 10 000 K. The las t colum n gives reference for th e adopted o scillator strength: 

NIST - iRalchenko et alj d2009r): K07 -iKurucz) d2007h : BBCB - iBerry et alj dl97lT> : S R - ISmith & Raggettl dl98lh : S -ISmithl dl988h; T - 
iTheodosioul d 19891) ; MFW - iMartin et alj dl988ri; BWL - lO'Brian et all d!991l) ; BKK - iBard et alj < 199 lb : HLGH - iHannaford et alj d!982h : 
BK- lBard& KocU dl994h : MW - iMiles & Wiesd dl969h . 
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Table 4. Abundances of individual elements. 



Ion 


Mines 


1 ,»,w AT / AT \ 

10gyVel/«totalJ 


Error 


1^.-./ AT 1 AT \ 

lOgfVel/'VtotaljQ 


Ci 


4 


-3.91 


±0.20 


-3.65 


Oi 


2 


-3.42 


±0.36 


-3.38 


Nai 


4 


-6.03 


±0.22 


-5.87 


Mg ii 


1 


-4.41 




-4.51 


Al i 


1 


-6.69 




-5.67 


S i 


1 


-4.76 




-4.90 


Sc ii 


6 


-8.97 


±0.59 


-8.87 


Ti i 


6 


-7.38 


±0.31 


-7.14 


Ti ii 


12 


-7.32 


±0.30 


-7.14 


V i 


8 


-7.15 


±0.41 


-8.04 


V ii 


6 


-7.28 


±0.32 


-8.04 


Cri 


12 


-6.45 


±0.34 


-6.40 


Cr ii 


16 


-5.91 


±0.30 


-6.40 


Mm 


8 


-6.26 


±0.29 


-6.65 


Mn ii 


3 


-5.68 


±0.20 


-6.65 


Co i 


40 


-6.00 


±0.22 


-7.12 


Con 


2 


-5.09 


±0.47 


-7.12 


Ni i 


8 


-7.01 


±0.47 


-5.81 


Cu i 


2 


-7.97 


±0.05 


-7.83 


Sri 


1 


-7.60 




-9.12 


Si" n 


2 


-7.91 


±0.07 


-9.12 


Y i 


1 


-7.97 




-9.83 


Y ii 


8 


-9.17 


±0.24 


-9.83 


Zri 


5 


-7.90 


±0.46 


-9.45 


Zr ii 


1 1 


-7.53 


±0.21 


-9.45 


Nb i 


4 


-7.48 


±0.25 


-10.62 


Mo i 


5 


-7.56 


±0.29 


-10.12 


Ru i 


4 


-7.57 


±0.34 


-10.20 


Rh i 


2 


-7.16 


±0.61 


-10.92 


Pd i 


6 


-6.98 


±0.23 


-10.38 


Cdi 


1 


-7.94 




-10.27 


So ii 


1 


-7.48 




-10.04 


La ii 


31 


-7.91 


±0.36 


-10.91 


Cen 


1 18 


-7.53 


±0.32 


-10.34 


Ce in 


4 


-5.70 


±0.24 


-10.34 


Pr ii 


37 


-8.54 


±0.31 


-11.46 


Prm 


26 


-6.58 


±0.38 


-11.46 


Ndi 


4 


-6.88 


±0.12 


-10.59 


Ndn 


109 


-7.1 1 


±0.34 


-10.59 


Ndm 


8 


-5.97 


±0.60 


-10.59 


Sm i 


1 


-7.1 1 




- 1 1 .04 


Smn 


63 


-7.25 


±0.39 


- 1 1 .04 


Eun 


7 


-8.08 


±0.45 


-11.52 


Gdn 


59 


-7.34 


±0.33 


-10.92 


Tbn 


3 


-8.84 


±0.22 


-1 1.76 


Tbin 


10 


-6.07 


±0.43 


-1 1.76 


Dy i 


1 


-7.88 




-10.90 


Dy ii 


29 


-7.45 


±0.33 


-10.90 


Dy in 


9 


-5.78 


±0.33 


-10.90 


Hoi 


1 


-6.79 




-11.53 


Ho hi 


4 


-6.34 


±0.11 


-1 1.53 


Er i 


1 


-8.16 




-11.11 


El ii 


23 


-7.81 


±0.30 


-11.11 


Erm 


4 


-6.18 


±0.29 


-11.11 


Tmn 


18 


-8.16 


±0.32 


-12.04 


Ybi 


1 


-7.44 




-10.96 


Ybii 


9 


-8.98 


±0.37 


-10.96 


Lu ii 


7 


-8.44 


±0.30 


-11.98 


Hf i 


1 


-7.89 




-11.16 


Hf ii 


3 


-8.39 


±0.30 


-11.16 


Ta i 


3 


-8.59 


±0.49 


-12.21 


Ta ii 


1 


-8.32 




-12.21 


V\ 1 


4 


—7 83 


±0 61 


in oi 


Wn 


3 


-7.79 


±0.35 


-10.93 


Re i 


2 


-7.60 


±0.04 


-11.81 


Ir i 


2 


-7.17 


±0.11 


-10.66 


Pti 


2 


-7.43 


±0.12 


-10.40 


Aui 


1 


-8.33 


±0.00 


-11.03 


Hgi 


2 


-7.76 


±0.15 


-10.91 


Thn 


9 


-9.18 


±0.35 


-11.98 


Thin 


5 


-8.09 


±0.18 


-11.98 


On 


17 


-8.44 


±0.15 


-12.56 



The solar abundances are taken from Gre vesse et 3] ((2007). 
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Fig. 8. Observed and theoretical energy distribution of HD 101065 calculated with scaled REE opacity from ISAN line list (see 
text). For comparison purpose the theoretical energy distribution of a normal F-type star HD 49933 is also shown. All theoretical 
calculations are smoothed with 20A Gaussian. Theoretical fluxes were scaled assuming radiuses: R = 1.98/? Q for HD 49933 and 
r eff = 6400 K models; R = 1 .90R Q for T eff = 6500 K; R = 1 .82/? for T eft - = 6600 K. 



